from math import *
#x=[24.533,16.033,35.95,7.85,0.583,16.867,6.95,0.333,18.25,10.583,14.8,7.333]
#y=[0.34,0.264,0.443,0.19,0.126,0.271,0.182,0.123,0.283,0.215,0.252,0.185]
print("régression lineaire y = a + b * x\n".center(40))
ch=input("entrez les valeurs x : ")
x=list(eval(ch))
ch=input("entrez les valeurs y : ")
y=list(eval(ch))
n=len(x)
sx=sum(x)
sx2=sum([w*w for w in x])
sy=sum(y)
sy2=sum([w*w for w in y])
sxy=sum([w*z for w,z in zip(x,y)])
a=(sy*sx2-sx*sxy)/(n*sx2-sx**2)
b=(n*sxy-sx*sy)/(n*sx2-sx**2)
d=(n*sxy-sx*sy)**2/((n*sx2-sx**2)*(n*sy2-sy**2))
print("{:s}{:>10.3g}{:s}{:>10.3g}{:s}{:>10.3g}{:s}".format("\na : ",
    a,"\nb : ",b,"\nd : ",d,"\n\n      x       y-calc   y-donné"))
for i in range(0,len(x)):
   
    print("{:>10.3f}{:>10.3f}{:>10.3f}".format(x[i],a+b*x[i],y[i]))
print()